AD-A033  033 


AERONAUTICAL  RESEARCH  ASSOCIATES  OF  PRINCETON  INC  N J F/6  20/11 

THE  SENSITIVITY  OF  IMPACT  RESPONSE  TO  THE  MECHANICAL  BEHAVIOR  O--ETC(U) 


I 


JAN  76  T B MCDONOUGH’  R M CONTILIANO  F44620-74-C-0045 

UNCLASSIFIED  ARAP-266  AFOSR-TR-76-1083  NL 


1 °F  1 

AD 

AC33033 

& 

I 

- _• 

1#  l 

/ 

i'll  ] 

u !;  Vi 

/®  i 

' 1; 

i 

] 

s 

; \(j 

1 A 

uT 

i'  All 

' 

- 

c! 

* / 

EI^C 

DATE 

FILMED 

2-' 

r 

n 

MICROCOPY  RESOLUTION  TEST  CHAflT 

NATIONAL  BUREAU  OE  STANDARDS-1963-j< 

' ^ f 

~~~  C-f 

-■  til 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wh»n  Data  Entaracf) 


READ  INSTRUCTIONS 


ICUMENTATION  PAGE 


BEFORE  COMPLETING  FORM 


REPORT NUMBEI 


S.  RECIPIENT'S  CATALOG  NUMBER 


GOVT  ACCESSION  NO 


LE  (and  Subtitle ) 


TVPF  OF  BCBMT  A 

Interim  rep 


VERED 


.THE  SENSITIVITY  OF  JMPACT  RESPONSE  TO  THE 
MECHANICAL  BEHAVIOR  OF  THE  PENETRATOR^-  ■< 


A.R.A.P.  Report  No.  266 


7.  author^; 


8.  CONTRACT  OR  GRANT  NUMBERfa) 


Thomas  B./ McDonough 
Ross  M. /Contiliand 


F4462)6-74-C 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 


9.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 


AREA  4 WORK  UNIT  NUMBERS 

681307 

9782-05 


Aeronautical  Research  Associates 

of  Princeton,  Inc.  / 

50  Washington  Road,  Princeton,  NJ  08540 


61102F 


<1.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


Air  Force  Office  of  Scientific  Research/NA" 
Building  410,  Bolling  Air  Force  Base 
Washington,  DC  20332 

u.  MONITORING  AGENCY  name  A AODRESSfif  different  from  Controlling  Office) 


15.  SECURITY  CLA! 


a report) 


Unclassified 


15a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  ti 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  abatract  entered  In  Block  20,  If  different  from  Report) 


18.  supplementary  NOTES 


19.  KEY  WORDS  (Continue  on  teverae  aide  It  neceeaary  and  Identify  by  block  number) 

Impact 

Integral  Theory 
Hydrodynamic 
Rigid-Plastic  Solid 
Penetration 


20.  ABSTRACT  (Continue  on  reverae  aide  It  neceeaary  and  Identify  by  block  number) 


An  integral  theory  is  under  development  at  A.R.A.P.  It 

Ts'*§  simple  theory/which  makes  It  easy  to  focus  attention  on  the 
various  physical  phenomena  which  occur,  and-/Tts*is  simple  enough  so 
that  the  performance  of  parametric  studies  is  economical?  The 
theory  has  progressed  In  steps,  starting  from  the  simplest  model 
possible  and  gradually  adding  more  complexities.  This  report 
presents  the  results  of  one  step  in  the  development:  a study  of 


EDITION  OF  I NOV  «s  1$  obsolete 


UNCLASSIFIED  O C><5 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fOTiwt  Dmla  Entered) 


SECURITY  CLASSIFICATION  OF  THIS  PAOE(TTh«n  Pti  l»HnO 


20.  (cont'd) 

? the  sensitivity  of  impact  characteristics  to  the  properties  of 
the  impacting  particle.  Four  ideal  models  of  behavior  are 
treated:  hydrodynamic,  Newtonian  fluid,  rigid-plastic  solid,  and 

rigid.  The  effects  on  the  impact  process  are  presented  in  terms 
of  the  partitioning  of  the  energy  between  particle  and  target, 
the  deformation  of  the  impacting  particle,  depth  into  the  target 
the  particle  penetrates,  and  the  target  mass  lost  as  a result  of 
the  cavity  formation.  The  predicted  trends  with  increasing 
impact  velocity  show  some  interesting  and  unforeseen  results; 
for  example,  penetration  depth  does  not  increase  always  with 
increase  of  particle  strength. 

\ 


NOMENCLATURE 


D 


b 
C 
e 
E 
K 
1 

*o 

mp 

mt 

R 

s 

t 

v± 

Vcm 

V* 

y 

Z, 


'ij 


ij 


4> 

n 


u 

Y 

(•) 


particle  width 

drag  coefficient 

deviator  of  stretching  tensor 

internal  energy  of  particle 

kinetic  energy  of  particle 

particle  height 

initial  particle  height,  width,  or  length 
particle  mass 

target  mass  removed  by  the  impact 

radius  of  curvature  of  the  flowing  shocked  layer 

axial  component  of  stress  deviator 

time 

speed  of  particle  front 

speed  of  the  particle's  mass  center 

characteristic  velocity  of  the  target 

penetration  depth 

particle  velocity 

stretching  tensor 

density  of  particle 

density  of  target 

stress  tensor 

strength  of  target  normal  to  surface 
strength  of  target  parallel  to  surface 
yield  strength  of  particle 

incidence  angle  along  surface  of  shocked  layer 

rectangular  coordinates  originating  at  particle's  mass 
center 

viscosity  of  liquid  droplet 
surface  tension 


iii 


I . INTRODUCTION  1 

II.  DEVELOPMENT  OF  INTEGRAL  THEORY  2 

Particle  Motion  2 

Cavity  Formation  7 

Dimensionless  Equations  12 

III.  RESULTS  OF  PARAMETRIC  STUDY  15 

A.  Energy  Partitioning  16 

B.  Particle  Shape  45 

C.  Depth  of  Penetration  52 

D.  Target  Erosion  61 

IV.  CONCLUSIONS  72 

REFERENCES  73 


1 


i 


I.  INTRODUCTION 


There  is  considerable  current  interest  in  high-speed  impact 
phenomena,  and  numerous  studies  presently  underway  are  being 
reported.  Generally,  these  studies  follow  one  of  two  approaches: 
on  the  one  hand,  experimental  programs  are  conducted  and  data  are 
extrapolated  empirically;  on  the  other  hand,  theoretical  studies 
have  led  to  the  development  and  use  of  very  detailed  multi-element 
computational  codes.  A third  approach  - the  development  of  an 
approximate  theory  - has  received  only  minimal  attention. 

The  development  of  a simple  theory  has  been  initiated  at 
A.R.A.P.  The  goal  of  this  theory  is  to  identify  the  way  in  which 
the  physical  properties  of  the  target  and  penetrator  are  related 
to  the  impact  phenomena.  Because  it  is  a simple  theory,  this 
approach  makes  it  easy  to  focus  one's  attention  on  the  importance 
of  any  one  of  the  many  physical  phenomena  that  enter  the  problem, 
such  as  the  energy  partition  between  particle  and  target.  Numerical 
computations  are  economical.  Thus,  parametric  studies  can  be 
conducted  reasonably  and  are  useful  for  interpreting  test  observa- 
tions in  terms  of  the  physical  processes  at  work  and  for  guiding 
the  design  of  effective  armor  systems. 

The  integral  theory  bears  a similar  relationship  to  the  multi- 
element computational  codes  that  the  integral  theory  of  boundary 
layers  bears  to  finite  difference  solutions.  It  Is  a simplified 
theory  based  on  heuristic  procedure.  The  development  has  followed 
a step-by-step  process,  starting  from  a very  simple  form,  and 
further  complexities  have  been  added  only  as  required  to  bring  the 
predictions  into  better  agreement  with  experimental  observations. 

One  aspect  of  the  impact  phenomenon  which  has  been  addressed  is 
characterizing  the  state  of  the  impacting  body.  Previous  studies 
have  shown  surprisingly  good  correlation  with  available  test  data1 
for  two  extreme  representations  of  penetrator  behavior  - rigid  and 
hydrodynamic.  In  the  present  study,  other  models  of  behavior  are 
introduced  which  include  material  strength  and  viscosity.  Finally, 
the  results  of  a parametric  study  are  presented  which  illustrate  the 
sensitivity  of  pertinent  impact  parameters  to  the  behavior  of  the 
impactor . 
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II.  DEVELOPMENT  OF  INTEGRAL  THEORY 


Particle  Motion 

We  start  with  the  global  conservation  equations  for  linear 
momentum  and  kinetic  energy  of  a fixed  body,  i.e.,  a material 
volume 


r unj da  dv  ■ It  foi i 

•'S  J\r  J \r 


Ji i°ijVa  + | dv  = £t  f 


zizi  d 

p 2 dv  + dt  E 


: represents  the  internal  energy  of  the  body,  and  its  prescrip- 

tion limits  the  energy  transfer  permitted  in  the  theory.  In  the 
present  case,  we  assume  the  impact  occurs  adiabatically , and  we 
characterize  the  energy  transfer  totally  in  terms  of  the  stress 
power  and  free  surface  energy.  The  latter  is  represented  only  in 
terms  of  the  equilibrium  properties.  Therefore, 


dt  E ~ J °ijZi,t 

v v 


dv  + I y da 


The  above  equations  are  complemented  by  the  conservation  of  mass 


a_  r 

dt  I 

J -n 


p dv  = 0 


and  a constitutive  law.  Because  elastic  strains  are  so  very  small 
compared  to  strain  associated  with  flow,  they  are  neglected.  The 
constitutive  law  is  represented  by  four  ideal  models  of  behavior: 

Hydrodynamic  model 

S1J  ' 0 (5! 
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Newtonian  fluid  model 


su  ■ 2lieij 


(5b) 


Rigid  plastic  model 


SiJ  SiJ 


4 2 
3 ay 


for  e^j  / 0 


(50 


Rigid  model 


£u  ' 0 


(5b) 


The  above  equations  represent  the  physics  of  particle  behavior 
in  general  form.  Further  development  of  the  present  theory  intro- 
duces no  new  physics.  Rather,  assumptions  are  introduced  into  the 
above  equations  which  reduce  them  to  an  approximate  set  of  ordinary 
differential  equations,  in  the  spirit  of  integral  boundary  layer 
theory  and  finite  element  analysis.  In  all  of  these  approximations, 
we  are  careful  to  satisfy  the  global  conservations  laws  for  the 
particle . 

Now  consider  the  normal  Impact  of  the  particle  onto  a semi- 
infinite target.  We  characterize  the  particle  as  a rectangular 
parallelepiped  of  square  planform.  For  simplicity,  the  deformation 
is  constrained  so  that  the  particle  remains  a rectangular  parallel- 
epiped. Penetration  into  a target  is  schematically  Illustrated  by 
Figure  1.  The  only  force  acting  on  the  particle  Is  the  contact 
force  with  the  target.  Figure  2 illustrates  a "free  body"  diagram 
of  the  particle  with  the  resistive  force  acting  to  Impede  penetra- 
tion. Regarding  the  flow  of  material  within  the  particle,  we 
assume  that  the  velocity  varies  linearly  with  position,  and  we 
neglect  compressibility.  The  situation  is  Illustrated  by  Figure  3- 
Coordinates  £ , n , x,  are  measured  from  the  mass  center,  and  the 

material  velocity  is  denoted  by  V , V , V . The  velocity  field 

x y z 

is  written  simply  in  terms  of  the  front  face  velocity  (Vj  and  the 

mass  center  velocity  (V  ): 

cm 


x \ b/2 / 2 
Vy  ■ Vcn,  +(l7?)( 
VZ  - (bT?)?* 


V,  - V 
1 cm , 


V,  - V = i-  Jl 
-L  cm  2 


9V  9 V 9 V 

— — + —X  + — - = 0 

9x  9y  9z  u 

The  conservation  equations  can  be  simplified  by  substituting 

(6)  and  (3)  into  (1)  and  (2): 


' DRAG  = Vcm 


DRAG  • V.  = K + E 


m P o 

v2 

+ l fb2  + H2  ' 

2 cm 

6 \ 2 

o n • 
f b si 

+ 2y  (i  - r)bi 

2 „ 

p b X. 

P 

and  substituting  (6)  and  (7)  into  (8)  gives  the  following: 


I ■ 


The  system  of  equations  composed  of  (7),  (9),  (10),  (11),  and 
one  of  the  models  represented  by  (5)  is  a closed  system  representing 
the  equations  governing  the  motion  of  the  particle.  The  drag  force 
represents  the  resistance  the  target  offers  to  penetration.  As 
previously,  we  assume  the  particle  penetrates  as  long  as  the 
Newtonian  pressure  at  the  front  face  is  greater  than  the  target 
strength  (called  o^).  While  penetration  continues,  the  drag  is 
assumed  proportional  to  the  Newtonian  pressure: 


if 


DRAG  = 


ptvi  > 


W 


„2W2 

PtvAb 


(12) 

(13) 


Consistent  with  (13),  a characteristic  velocity  can  be  defined  for  a 
target,  below  which  penetration  ceases. 


(14) 


V = 

The  equations  for  penetration  are  now  complete.  Next  we  address 
mass  removal. 


Cavity  Formation 

We  visualize  the  flow  of  target  material  around  the  particle 
as  illustrated  by  Figure  4.  There  is  a disturbed  layer  of  target 
material  in  front  of  the  particle.  Fixing  ourselvqs  to  the  front 
face,  we  see  the  target  approach  at  velocity  VA  (see  the  figure), 
and  we  assume  the  material  flows  laterally,  leaving  the  particle 
edge  at  the  same  velocity  VA  . Neglecting  compressibility,  the 
mass  flowing  into  balances  that  flowing  out  of  the  disturbed  region 
if  the  thickness  is  b/4  as  illustrated. 

The  two  layers  (particle  and  disturbed  target)  possess  signi- 
ficant lateral  momentum  in  their  flow  outward.  It  Is  this  momentum 
that  cuts  a wide  cavity  in  the  target.  The  flow  does  not  continue 
flat;  it  Is  turned  upward  by  the  target  pressure  acting  on  it  from 
below.  At  a particular  section,  we  can  estimate  the  curvature. 
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Figure  4.  Model  of  cavity  formation 


Consider  the  momentum  delivered  to  the  disturbed  layer;  the 
pressure  equals  the  net  moment. 


Ap  = p cos2  <J>  ( 

But  this  pressure  is  what  produces  the  rate  of  change  of  momentum: 


2 2 

p VA  cos  <|>  = 


If  we  assume  the  lateral  velocity  is  constant  through  the  thickness 
of  each  layer  independently,  the  integral  has  the  following  approxi- 
mate solution: 


J dz  eP(V2  *n  fwT77!d+  ptvt  *n 


When  this  result  is  used  in  (16),  the  following  solution  can  be 
derived : 


L\OOS*J  'UAvJ 


4R  i \ / HR 

4R  V~!nran  brr 

b ' 1 ' b / \b 


We  use  the  value  at  the  particle  edge  for  vj/vt  » and 


everywhere 


because  both  V and  V,  vary  linearly  with  distance  from  the 
P ^ 

center.  Finally,  the  equation  for  the  radius  of  curvature  becomes 


2.  M b/2 

cos  ♦ = Ur  hr 


hr  _ 1 \ Hr 

HR  b 57  + ( 4Rb  , 

TT  “ b / \ b _ 1 


To  consider  the  flow  of  disturbed  target  layer  outside  of  the 
particle  edge,  one  needs  to  neglect  the  lateral  momentum  of  the 
particle  in  the  above  equations. 


The  lateral  momentum  of  the  flow  will  continue  to  cut  a wider 
cavity  until  the  transverse  strength  of  the  target  is  great  enough 
to  stop  the  flow.  To  analyze  this  situation,  we  assume  the  flow 
will  be  turned  up  at  an  angle  6 by  a rigid  wall  and  illustrate  it 
by  Figure  5.  The  Newtonian  pressure  in  the  flowing  layer 
(produced  by  the  momentum  flux  from  the  undisturbed  region)  at  the 
angle  0 is 

Pq  = p^V^  cos20  (2; 

To  turn  the  flow  vertical  requires  a pressure  Increase  equal  to  the 
horizontal  momentum  flux  of  both  layers  : 


Jlp  (b/2)2  + (b/4)p  V2  ? 

AP  = — 2 k cos  0 

U + b/4 ) 


Therefore,  the  wall  must  sustain  a pressure  of 


pl  = Po  + AP  = ! 


/ \ ( ^ • 2 

(*  * *177 ) ptv^  + (T7TT)  COi 


or  it  cannot  turn  the  flow.  If  we  equate  this  pressure  to  the 
transverse  strength  of  the  target  (o2),  we  can  solve  for  0^^  , the 
angle  at  which  the  flow  will  be  turned 

-i  r °2  i1/c 

ei  = cos  7-T7  (25) 
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where 


f = 


1 + 


-77-  + 1 

b 


(£) 


»tvi 


The  shape  of  the  cavity  that  is  created  by  the  flow  at  a 
typical  instant  during  the  penetration  is  illustrated  by  Figure 
6a.  The  final  shape  of  the  cavity  is  shown  by  Figure  6b.  As 
suggested  by  the  figure,  it  is  determined  by  superposing  all  the 
shapes  computed  at  many  time  steps  during  the  penetration  and 
using  their  envelope. 


Dimensionless  Equations 

It  Is  informative  to  define  dimensionless  parameters.  In  fact, 
that  is  an  advantage  of  an  approximate  theory.  The  problem  is 
reduced  to  a manageable  form  so  that  "key"  dimensionless  parameters 
can  be  identified.  In  this  report,  barred  quantities  will  denote 
dimensionless  parameters  whose  definitions  are  given  in  the  nomen- 
clature. When  the  governing  equations  are  cast  in  dimensionless 
terms,  they  become  the  following: 


Penetration 


v'  = - y PV2b2 
cm  2 

K'  + E’  = - CDp  v]b2 


V2  + ~ (n2 
cm  6 


2l2) 


(26) 


V.  - V 


= 6pb2s£  + 4 (l  - — ) byC 
' l ' 

= l 


cm 

b2I  = 1 


= ! V 
- 2 * 


n = b ' 
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Hydrodynamic 


Constitutive  Law 


4p  - 
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Newtonian  Fluid 
Rigid  Plastic 
Rigid 


Cavity  Formation 


1 

/db  Y 

1 

\dx  j 

„2 

p 

4R  _ 1 _ 


+ £n  -= 

1 


2 

cos  4> 


= Hn 


az- 1 


where 
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III.  RESULTS  OP  PARAMETRIC  STUDY 


The  results  of  a parametric  Investigation  of  impact  cratering 
by  high-speed  particles  are  presented  in  this  chapter.  The  results 
describe  one  particular  class  of  impacts,  namely,  normal  impact  of 
a single  particle.  An  understanding  of  the  impact  process  for  this 
class  is  a prime  requisite  for  the  more  usual  but  more  complicated 
classes  which  include  multiple  and  oblique  impacts.  Hence,  the 
results  which  follow,  aside  from  their  obvious  usefulness  for  data 
interpretation,  parameter  sensitivity  studies,  and  test  planning, 
will  be  used  to  illustrate  many  of  the  basic  physical  mechanisms 
which  control  the  impact  process. 

The  impact  conditions  which  were  used  for  this  study  encompass 

a wide  range  of  particle  and  target  density  and  impact  velocity. 

Three  values  of  the  ratio  of  target  density  to  particle  density  (p) 

were  considered:  p = 0.1  (i.e.,  a "heavy"  particle),  p = 1.0  , 

and  p = 10.0  (i.e.,  a "light"  particle).  For  each  value  of  p , 

— * 
the  value  of  V (the  ratio  of  particle  Impact  velocity  to  V , 

the  velocity  characteristic  of  the  target  strength)  was  varied  from 

slightly  above  1.0  to  approximately  50.  In  the  integral  theory, 

V = 1.0  represents  the  velocity  ratio  for  incipient  target  damage; 

Impact  velocities  below  this  value  produce  elastic  response  in  the 

target  but  produce  no  penetration  and  are,  therefore,  excluded  from 

this  analysis. 

Two  mechanisms  are  included  in  the  integral  theory  to  account 
for  energy  dissipation  within  the  particle,  viscous  dissipation  for 
liquid  particles  and  plastic  dissipation  for  solid  particles.  Para- 
metric variation  of  the  Internal  energy  dissipation  was  achieved 
using  particle  viscosity  (or  yield  strength)  values  from  zero 
(I.e.,  the  hydrodynamic  limit)  to  infinity  (I.e.,  the  rigid  limit). 

Many  parameters  are  needed  to  describe  the  impact  process; 
four  of  these  are  given  detailed  consideration  below.  They  include 
energy  partitioning  within  the  particle,  particle  deformation  and 
depth  of  penetration,  and  target  mass  removal. 


A.  Energy  Partitioning 

Prior  to  impact,  the  total  energy  contained  in  the  system 

comprising  the  particle  and  target  is  the  kinetic  energy  of  the 

particle*  (K  ).  During  an  impact,  work  is  done  on  the  target  and 

energy  is  transferred  at  the  expense  of  K . In  this  study,  four 

modes  are  accounted  for:  kinetic  energy  of  the  particle  associated 

with  mass  center  motion  (K^  ),  kinetic  energy  of  the  particle  assoc- 

cg  — 

iated  with  motion  relative  to  the  mass  center  (KR),  energy  that  is 
dissipated  by  the  resistance  of  the  particle  to  flow  (E),  and  work 
done  on  the  target  (E^).  Although  the  energy  may  be  transferred 
from  one  mode  to  another,  it  must  remain  In  the  system  somewhere. 

Its  conservation  requires  that  the  energy  components  obey  the 
following  equation  at  any  time: 

K = K + K0  + E_  + E (28) 

o eg  R D 

A study  of  how  the  energy  partitions  for  various  impact  cases  is 
very  enlightening  for  the  impact  process. 

Rigid  Particle 

The  rigid  particle  idealization  is  also  approached  in  the  limit 
by  very  strong  solid  particles  and  very  viscous  liquid  particles  for 
an  appreciable  range  of  impact  velocity.  By  definition,  the  rigid 
particle  does  not  deform.  Hence  it  has  no  relative  kinetic  energy 
and  dissipates  no  energy  internally.  For  this  case,  the  energy 
equation  becomes 

K = K + En  (rigid  particle)  (29) 

o eg  u 

Since  the  velocity  is  uniform  within  the  particle,  it  follows  that 


,1 


- K V2  V2 

JL.  = -££.  = .eg.  = _A 

K K V 2 V2 

o o o o 

I.e.,  the  fraction  of  the  initial  kinetic  energy  which  remains  in 


(30) 


*The  internal  energy  of  both  the  particle  and  target  is  taken  to  be 
zero  prior  to  impact.  Internal  energies  are  then  measured  relative 
to  this  ground  state. 
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the  particle  at  any  instant  during  impact  is  simply  the  square  of 
the  ratio  of  particle  velocity  to  impact  velocity.  In  particular, 
at  the  end  of  impact,  VA  = 1 and,  therefore. 


- K 
K _ eg 


K 


K 


o o 

Equation  (28)  requires  that 


1 

n 


(end  of  impact) 


(31) 


i . I 

K V2 

o o 


(end  of  impact) 


(32) 


Note  that  the  energy  partitioning  for  a rigid  particle  is  independent 
of  the  density  ratio. 

Figure  7 shows  the  energy  partitioning  for  a rigid  particle 
according  to  Eqs.  (3D  and  (32).  This  figure  and  all  subsequent  ones 
of  this  type  show  the  energy  partition  at  the  end  of  Impact  as  a 

function  of  the  initial  Impact  velocity.  The  end  of  impact  is 

* _ 

defined  as  the  instant  when  V,  = V (i.e.,V,  =1).  A method  for 

* 

estimating  the  characteristic  velocity  V is  described  in  Refer- 
ence 1.  Although  the  particle  still  has  axial  velocity  at  this 
Instant,  the  kinetic  energy  associated  with  further  deceleration  of 
the  particle  is  stored  elastically  in  the  target  and  causes  no 
additional  penetration  of  the  target. 

Figure  7 shows  a very  rapid  increase  in  the  transfer  of  energy 
from  that  stored  in  elastic  response  to  work  done  on  the  target  as 
the  impact  velocity  is  increased.  The  latter  component  dominates 
for  velocity  ratios  above  1.5.  Note  that  nearly  all  of  the  particle 
kinetic  energy  is  used  to  damage  the  target  when  the  velocity  ratio 
is  larger  than  10. 


Hydrodynamic  Particle 

The  hydrodynamic  particle  is  another  useful  idealization.  By 
definition,  it  is  a particle  which  deforms  but  which  offers  no 
resistance  to  shearing  flow.  It  is  a fairly  good  approximation  for 
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Figure  7.  Energy  partition  vs.  impact  velocity 


a low  viscosity  liquid  such  as  water.  The  hydrodynamic  particle 
introduces  a third  component  into  the  energy  partition,  i.e., 
the  relative  kinetic  energy  due  to  deformation.  Thus, 


K = K + K„  + En  (hydrodynamic)  (33) 

o eg  n V 


Figures  8,  9,  and  10  show  the  energy  components  as  a function 
of  impact  velocity  for  density  ratios  of  0.1,  1.0,  and  10, 
respectively.  Several  trends  should  be  noted.  First,  for  low 
impact  velocities,  most  of  the  energy  is  stored  as  elastic  response 
of  the  target.  The  lighter  the  particle,  the  larger  the  velocity 
at  which  this  component  dominates. 

Second,  for  very  large  velocities,  the  elastic  energy  is 
insignificant  compared  to  both  the  lateral  kinetic  energy  of  the 
particle  and  the  work  performed  on  the  target.  For  the  heavy 
particle,  90$  of  the  energy  represents  work  done  on  the  target  and 
10$  represents  the  lateral  kinetic  energy  in  the  particle.  In 
contrast,  the  light  particle  deforms  more  and  does  less  work  on  the 
target.  As  a result,  the  two  components  are  nearly  equal. 

Third,  another  difference  between  high  and  low  p Impacts 
should  be  noted  for  velocity  ratios  above  10.  When  p < 1 , the 
particle  is  fully  embedded  in  the  target  (y/I  > 1).  Hence,  all  of 
the  energy  is  expended  on  the  target.  In  contrast,  when  p > 1 , 
the  particle  is  only  partially  embedded  in  the  target,  and  much  of 
the  energy  accompanies  the  particle  ejecta  and  is  dissipated  by  the 
air  surrounding  the  target. 

Transfer  of  energy  between  the  various  modes  in  the  transition 
regime  (2  < V < 10)  is  apparently  a complicated  process  which  a 
study  of  Figures  8-10  indicates.  A parameter  which  is  useful  for 
an  understanding  of  this  coupling  is  the  deformation  parameter  i,  . 
It  is  defined  as 


"Id*, 
K 2 dr 


(34) 


where  l represents  the  nondimensional  thickness  of  the  particle. 
In  the  absence  of  any  strong  elastic  restoring  forces,  the  particle 
is  compressed  during  impact.  Hence,  £ is  a negative  quantity. 
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gure  9.  Energy  partition  vs.  impact  ve 


Hydrodynamic  particle 
p=  10.0 


In  the  following  discussion,  it  is  convenient  to  recast  £ in 
terms  of  particle  velocity  and  to  work  with  positive  numbers. 

It  is  easy  to  show  that  for  a linear  velocity  field  the  deformation 
parameter  may  be  written 

- 5 - Vog  - (35) 

i.e.,  it  represents  the  rate  at  which  the  particle  center  of  mass 
approaches  the  impact  surface. 

Figure  11  shows  the  deformation  parameter  as  a function  of 
time  for  p = 1.0  and  for  several  values  of  impact  velocity.  The 
final  point  on  each  curve  coincides  with  the  end  of  impact.  It  is 
evident  that  the  initial  deformation  rate  is  greater  for  greater 
impact  velocity.  However,  during  a given  Impact,  the  rate  of 
deformation  decreases;  the  particle,  in  effect,  becomes  "stlffer." 
Eventually,  if  the  impact  velocity  Is  high  enough,  the  deformation 
parameter  reaches  a peak  value  and  then  decays.  For  the  case 
p = 1.0,  the  peak  Is  reached  during  impact  for  V > 3-  This  value 
is  the  same  as  the  velocity  at  which  the  sudden  transition  occurs 
in  Figure  9.  A similar  study  of  deformation  histories  for  p = 10 

indicates  flattening  at  V.  = 6.5,  which  also  corresponds  to  the 

A — 

transition  on  Figure  10.  For  p = 0.1,  the  study  shows  the  deforma- 
tion history  occurs  even  for  the  lowest  impact  velocities  treated. 

Further  illustration  is  contained  in  Figure  12  which  shows  the 
deformation  parameter  for  each  density  ratio  at  V = 3.32  . Three 
characteristics  should  be  noted.  First,  the  time  scale  of  the  impact 
process  varies  inversely  as  p , i.e.,  the  heavier  the  particle,  the 
more  time  required  to  decelerate  it.  Second,  the  peak  value  for  (-£) 
has  not  been  reached  for  the  light  particle,  whereas  it  has  been 
passed  for  the  heavy  particle.  Hence,  the  relative  kinetic  energy 
is  still  increasing  rapidly  for  the  light  particle  (Fig.  10)  whereas 
it  is  approaching  the  high  velocity  limit  for  the  heavy  particle 
(Fig.  8).  Finally,  the  value  of  (-£)  at  the  end  of  impact  increases 
with  p . The  latter  effect  helps  to  explain  the  curves  in 

Figures  8,  9,  and  10.  It  is  easy  to  show  that,  using  Eq.  (35), 
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r ™ 


Kcg  = (Vj.  - O (36: 

Since  Vj.  = 1 at  the  end  of  impact,  it  follows  that  as  (-£) 

increases,  K must  increase.  Hence,  the  value  of  K (i.e.,  the 

eg  c6_ 

elastic  energy  stored  in  the  target)  must  increase  with  p and  must 
exhibit  a rapid  change  in  slope  when  (-£)  peaks. 


Viscous  and  Plastic  Particles 

For  the  general  case  of  a viscous  or  plastic  particle,  all  of 
the  components  in  Eq.(28)  must  be  considered.  Two  models  were 
considered  for  the  internal  energy  dissipation:  viscous  and  plastic 

deformation.  Figures  13  through  23  show  the  energy  partitioning  for 
a viscous  particle;  Figures  2A  through  29  show  corresponding  results 
for  a plastic  particle.  Each  figure  shows  the  energy  components  for 
a range  of  impact  velocities  and  for  a particular  combination  of  o 
and  y (or  a ) . 

Figures  13  through  16  are  for  p = 1 and  a range  of  y from 


0.01  to  10.  The  results  for  this  density  ratio  are  typical  of  the 

results  shown  for  the  other  density  ratios.  As  expected,  for  low 

viscosity  (Fig.  13)  the  energy  partitioning  is  similar  to  the 

hydrodynamic  results  of  Figure  9.  The  internal  energy  dissipation 

never  exceeds  3%  of  the  initial  particle  kinetic  energy. 

Figure  14  shows  the  results  for  a larger  value  of  viscosity 

(y  = 0.1).  For  V = 4 , almost  20%  of  the  particle  energy  is 

dissipated  as  heat.  Before  proceeding  further,  it  should  be  noted 

that  it  is  the  ratio  of  E/K  which  has  a maximum  value  as  a 

o 

function  of  impact  velocity.  The  energy  dissipation  itself  is  a 
monotonically  increasing  function  of  impact  velocity.  Hence,  E 
does  not  approach  zero  for  large  velocities.  On  the  contrary,  it 
is  growing  quite  large;  it  is  just  increasing  more  slowly  than  K. 
Both  the  maximum  value  for  E/Kq  and  the  velocity  at  which  the 


maximum  occurs  Increase  directly  with  y 


For  example,  when  y = 1 


(Fig.  15),  the  maximum  value  of  E/Kq  is  25%  and  it  occurs  when 
V = 35;and  when  y = 10  (Fig.  16),  the  maximum  has  not  been  reached 
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Lgure  13.  Energy  partition  vs.  impact  velocity 


Viscous  particle 


Viscous  particle 


Figure  17.  Energy  partition  vs.  impact  velocity 


Viscous  particle 


Energy  partition  vs.  impact  velocity 


Viscous  particle 


Figure  26.  Energy  partition  vs.  impact  velocity 


Plastic  particle 


Figure  28.  Energy  partition  vs.  impact  velocity 


Plastic  particle 
> = 10.0  = 20.0 


gurc  29-  Energy  partition  vs.  impact  velocity 


for  V = 50.  For  both  these  conditions  of  relatively  high  viscosity, 
there  is  an  insignificant  amount  of  relative  kinetic  energy  in  the 
particle;  the  particle  has  simply  not  deformed  very  much  due  to  the 
viscous  stresses.  These  stresses  not  only  generate  more  heat  as 
noted  above,  but  the  "sharper"  particle  is  able  to  penetrate  deeper 
into  the  target  and  do  more  work  on  it.  This  additional  work  is 
shown  by  the  overshoot  of  the  ED  curves  in  Figures  15  and  16 
relative  to  the  hydrodynamic  limits. 

Figures  13  through  16  show  that  for  sufficiently  high  veloci- 
ties, the  energy  partitioning  approaches  the  hydrodynamic  values  of 
Figure  9.  This  is  clearly  evident  in  Figures  13  and  14.  In  Figure 
15,  Ep  is  approaching  the  limiting  value  of  70?,  but  a transfer 

between  E and  K„  must  still  occur  for  larger  impact  velocities. 

n 

Still  larger  velocities  are  required  before  the  limiting  conditions 
are  reached  as  in  Figure  16. 

Figures  17  through  20  are  for  p = 0.1  and  Figures  21  through 
23  are  for  p = 10.  The  results  are  similar  qualitatively  to  those 
described  above  for  p = 1.0. 

Figures  24  through  29  show  the  energy  partitioning  in  'a 
plastic  particle  for  all  combinations  of  the  three  density  ratios 
and  two  values  of  particle  yield  strength.  These  results  are 
similar  to  the  viscous  results  with  one  exception.  A viscous 
particle  deforms  continuously  under  an  applied  stress;  a plastic 
particle  deforms  only  when  the  applied  stress  exceeds  the  yield 
strength  of  the  particle  material.  It  can  be  shown  that  the  impact 
velocity  for  incipient  particle  deformation  (V  ) is  related  to 
particle  yield  strength  by 


= »^2o 


( 37) 


For  Impact  velocities  below  Vy  , there  is  no  particle  deformation, 
and  the  energy  partitioning  is  that  of  a rigid  particle  which  was 
described  earlier  in  Figure  7. 

The  constraint  imposed  by  Eq.  (37)  is  evident  in  Figures  24 
through  29.  For  example,  o =20  for  Figure  25  and,  therefore, 
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Vy  = 6.3.  In  the  Interval  1.0  < V < 6.3,  Figure  25  is  a replica 
of  Figure  7.  Due  to  Its  rigidity,  the  particle  does  more  work  on 
the  target  than  a deforming  particle  as  evidenced  by  the  overshoot 
of  the  ED  curve.  As  the  Impact  velocity  increases  above  V , 
the  particle  yields.  It  does  less  work  on  the  target  as  larger 
percentages  of  the  energy  are  expended  laterally  and  dissipated 
internally.  Eventually,  for  high  enough  impact  velocities,  the 
particle  strength  is  insignificant,  and  the  energy  partitioning 
again  approaches  the  hydrodynamic  values. 

B.  Particle  Shape 

Figures  30,  31,  and  32  show  particle  width  at  the  end  of 
impact  as  a function  of  impact  velocity  and  particle  viscosity 
for  three  values  of  density  ratio.  The  most  striking  feature  of 
these  curves  is  the  (apparent)  transition  of  behavior  between  low 
and  high  velocity  impacts.  The  transition,  which  is  most  apparent 
for  p = 1 and  10  and  for  small  values  of  u , corresponds  to  the 
stiffening  of  the  particle  which  was  discussed  In  Section  III. A. 

It  is  also  interesting  that  the  sharpness  of  the  transition  is 
blunted  by  increasing  viscosity. 

With  the  exception  of  impact  velocities  below  the  transition 
velocity,  the  hydrodynamic  solution  provides  an  upper  limit  for 
the  lateral  growth  of  the  particle.  At  large  impact  velocities 
(depending  on  the  value  of  p),  the  width  of  each  particle 
approaches  the  hydrodynamic  limit  closely. 

Figures  30  through  32  show  that  for  a given  particle  density 
and  impact  velocity,  viscosity  has  a larger  effect  on  particle 
shape  as  the  target  density  increases.  For  example,  for  p = 1 
and  V = 10,  the  width  of  the  particle  is  within  15%  of  the  hydro- 
dynamic  value  when  p = 0.1  , but  the  width  is  less  than  half  the 
hydrodynamic  value  when  p = 10.  The  implications  of  the  result 
as  regards  mass  erosion  are  noted  later. 

Figures  33,  3^,  and  35  show  an  analogous  set  of  particle 
width  curves  for  plastic  particles.  These  curves  exhibit  the  same 
trends  as  do  Figures  30  through  32  for  viscous  particles  - with 
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Viscous  part 


Figure  32.  Particle  width  vs.  impact  velocity 
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one  exception.  The  plastic  particle  does  not  deform  <nt 
yield  strength  is  exceeded.  As  a result,  plastic  j r*  it.  • u 
remain  rigid  for  the  duration  of  impact  for  low  vel  citit:  . 

For  somewhat  higher  velocities,  the  impact  pressure  c 
exceed  the  yield  strength,  and  the  particle  begins  to  defer:: 
However,  the  particle  may  become  rigid  during  impact*  as  a 
of  the  stiffening  behavior  described  earlier.  The  behavior 
represented  by  Figures  33  through  35  is  different,  depend'.:  . 
the  impact  velocity.  There  are  three  regimes.  At  low  v< 
a particle  does  not  deform;  at  high  velocity  the  part  it 
continuously  throughout  the  penetration;  arid  there  is  an  i o’  • 
late  regime  for  which  the  particle  deforms  for  a time  a: 
rlgidlzes . 


C.  Depth  of  Penetration 

Figures  36  through  41  shew  the  maximur  depth  of  penet - 1 

multiplied  by  density  ratio  as  a function  of  impact  vel  oi 
particle  viscosity.  In  terms  of  these  parameters,  rl ;:a  t 
behavior  produces  a straight  line  on  semi-log  paper  according 
the  relation 


py 


max 


In  V 


For  this  analysis,  a value  of  C^/2  = 1 was  used.  The  r _ 
described  below  are  applicable  to  other  valuer,  of  the  drag 
cient  if  one  considers  the  ordinate  of  plots  tv  ce  C,.  2 py 
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it  is  interesting  that  penetration  depth  for  hydrodynamic  rl  • 
is  a lower  bound  to  the  other  predictions.  Much  of  the  in. 
behavior  is  similar  for  viscous  and  plastic  particles.  i.  .e 
the  three  velocity  regimes  are  apparent  in  the  plots.  At 
velocity  and  high  strength  (i.e.,  either  u or  o ),  pred:  \ • 

y 

penetration  depth  follows  (closely  for  viscous  behavi.  r md  ...  * 
for  plastic  behavior)  the  rigid  particle  predictions.  A1  f •• 


*The  criterion  for  rigid! zation  of  a flowing  particle  i. 
This  follows  the  requirement  that  internal  energy  pr  d . • 
positive,  which  is  violated  in  the  present  analysis  if 
from  zero  to  a negative  value. 
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Figure  36.  Penetration  vs.  impact  velocity 


Plastic  particle 


velocity,  predicted  penetration  depth  approaches  closely  the 
hydrodynamic  particle  predictions.  The  transition  from  the  rigid 
to  the  hydrodynamic  limits  occurs  at  intermediate  impact  velocity. 

The  nature  of  viscosity  is  such  that  the  abruptness  of  the  transi- 
tion exhibited  by  hydrodynamic  behavior  is  blunted  for  larger  u - 
Higher  for  plastic  behavior  tends  to  blunt  the  same  abruptness, 

but  the  nature  of  plasticity  also  introduces  an  abruptness  of  its 
own  into  the  graphs. 

A very  interesting  phenomenon  is  exposed  by  Figures  38  and  Al. 
Depth  of  penetration  for  p < 1 is  bounded  from  above  by  rigid 
particle  predictions.  But  this  is  not  the  situation  when  the 
particle  is  less  dense  than  the  target  (i.e.,  0 > 1).  Penetrations 
as  much  as  35%  greater  are  predicted,  and  the  impact  velocity  at 
which  the  peak  occurs  depends  on  the  value  of  p or  o^  . Also, 
a particle  of  greater  strength  may  penetrate  less  than  one  of 
lesser  strength.  Furthermore,  increasing  the  impact  velocity  of  a 
given  particle  design  may  result  in  lower  penetration.  All  of 
these  characteristics  are  somewhat  surprising  and  motivate  further 
discussion  of  the  underlying  phenomenon. 

In  the  context  of  the  present  theory,  the  depth  of  penetration 
is  the  time  integral  of  the  front  face  velocity.  So  if  VA  is 
plotted  as  a function  of  t for  a given  impact,  then  the  area  under 
the  curve  is  numerically  equal  to  the  depth  of  penetration.  Detailed 
studies  of  individual  cases  show  that  the  VA  history,  hence  also 
the  depth,  is  very  sensitive  to  the  particle  defoliation  history 
which,  in  turn,  is  controlled  by  strength. 

The  abruptness  of  the  transition  for  p = 1C  and  hydrodynamic 
behavior  may  be  viewed  further  in  terms  of  the  VA  - t traces. 

Traces  for  three  values  of  impact  velocity  are  shown  by  Figure  A2. 
Note  that  for  VA  < 6,  the  front  face  velocity  decreases  rapidly 
to  1 and  ceases  penetration.  But  at  Vx  = 7 , the  trace  becomes 
flat  at  a value  in  excess  of  1,  and  the  penetration  persists  a 
longer  time.  The  depth  of  penetration  (i.e.,  the  area  under  the 
trace)  takes  a sudden  jump  at  Vj_  = 7.  That  the  front  face  decel- 
eration ceases  or  even  reverses  seems  at  first  inconsistent  with 
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conservation  principles.  Several  cases  have  been  checked  and 
no  inconsistencies  were  found.  For  example,  Figure  43  demonstrates 
that  the  mass  center  velocity  decreases  monotonically  with  time 
even  though  the  front  face  may  speed  up. 

D.  Target  Erosion 

The  coupled  effects  of  energy  partitioning,  particle  deforma- 
tion, and  depth  of  penetration  yield  the  target-mass  erosion 
curves  shown  in  Figures  44,  45,  and  46  for  viscous  particles  and 
in  Figures  47,  48,  and  49  for  plastic  particles.  The  oidinate  of 
each  figure,  which  is  denoted  by  m , represents  the  ratio  of  target 
mass  loss  to  particle  mass.  A cursory  glance  at  these  figures 
reveals  the  strong  Influence  of  particle  viscosity  (or  yield 
strength)  on  target  erosion. 

Before  examining  these  results  in  detail,  it  should  be  noted 
that  the  details  of  the  erosion  process  are  the  weakest  link  in  the 
integral  theory  at  this  time.  In  part  this  is  due  to  a paucity  of 
the  type  of  detailed  data  needed  to  develop  a general  model.  But 
just  as  Important  Is  the  fact  that  most  of  the  early  effort  in  the 
development  of  this  theory  has  been  devoted  to  the  prediction  of 
penetration  ana  particle  deformation.  The  theory  has  now  reached 
the  stage  where  effort  toward  cavity  predictions  is  warranted,  and 
an  improved  theory  of  mass  loss  is  under  development  at  present. 

Mass  loss  predictions  using  the  present  theory  are  compared  to 
experimental  observations  to  gain  insight  into  how  to  improve  the 
theory . 

A comparison  of  Figures  44  through  49  shows  certain  common 
trends.  First,  erosion  varies  inversely  as  the  density  ratio; 
the  heavier  particle  produces  the  most  damage,  all  other  conditions 
unchanged.  Second,  erosion  increases  with  viscosity  (or  yield 
strength)  for  low  impact  velocities  although  the  amount  of  mass 
removed  is  extremely  small.  In  contrast,  erosion  varies  inversely 
as  viscosity  for  high  impact  velocities;  the  hydrodynamic  value  is 
the  upper  limit  for  target  erosion. 
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These  contrasting  results  reflect  the  complex  coupling  of 
the  particle  deformation  and  penetration  and  the  amount  of  energy 
contained  in  each  of  the  available  modes.  It  was  shown  earlier 
that  an  increase  in  viscosity  results  in  deeper  penetration  but  a 
narrower  particle  and,  for  some  velocities,  less  lateral  kinetic 
energy.  But  greater  lateral  kinetic  energy  produces  greater 
cavity  width.  It  follows  that  cavity  volume  is  a very  complex 
function  of  particle  strength  and  impact  velocity. 

The  complexity  of  the  phenomenon  is  Illustrated  well  by  the 
predicted  mass  loss  behavior  in  the  Intermediate  velocity  regime. 
For  example,  in  the  interval  6 <_  V <_  15  , Figure  48  shows  that  it 
is  possible  to  decrease  erosion  by  increasing  the  Impact  velocity 
when  = 20  . In  this  velocity  interval,  the  width  of  the 
particle  has  increased  by  80%  (see  Fig.  3*0  hut  the  depth  of  the 
cavity  has  decreased  by  nearly  25%  (see  Fig.  40).  In  addition,  the 
increase  in  particle  kinetic  energy  associated  with  higher  impact 
velocities  does  not  produce  proportionate  Increases  in  lateral 
kinetic  energy  of  the  particle  or  in  work  done  by  the  particle  on 
the  target.  Rather,  most  of  the  additional  energy  is  used  to  heat 
the  particle  (Fig.  25).  As  a result,  despite  the  fact  that  the 
width  of  the  particle  has  Increased,  the  width  of  the  cavity  has 
actually  decreased.  Figure  50  shows  the  final  cavity  shape  for  the 
two  velocities  which  bound  the  regime.  The  reason  for  the  dip  in 
the  erosion  curve  is  evident  from  this  figure. 

As  the  velocity  is  increased  beyond  tils  regime,  more  of  the 
kinetic  energy  shows  up  as  relative  kinetic  energy  and,  as  a 
result,  the  cavity  becomes  much  wider.  Her.ce,  the  erosion  curve 
again  Increases  with  further  increase  of  impact  velocity. 

In  summation,  perhaps  the  most  important  point  of  these 
results  is  the  recognition  of  the  uniqueness  of  each  impact  condi- 
tion. Blanket  statements  such  as,  "The  stronger  the  particle,  the 
larger  the  erosion,"  or  "Erosion  increase''  ■ 1th  impact  velocity" 
simply  are  not  always  true.  Data  extrapolation  or  test  planning 
based  on  such  statements  may  be  faulty  and,  when  the  ultimate 


objective  is  hardware  design,  even  dangerous.  Hopefully,  analyses 
such  as  this  will  be  able  to  bound  the  problem  and  eliminate  the 
need  for  all-encompassing  but  unproven  general!  zations . 


IV.  CONCLUSIONS 


The  first  steps  in  the  creation  of  an  integral  theory  for 
impact  cratering  by  high-speed  particles  have  been  presented  and 
used  to  perform  a parametric  investigation  of  single  particle, 
normal  impact.  Density  ratio,  impact  velocity,  and  dissipation 
parameter  (particle  viscosity  or  yield  strength)  were  varied  over 
a wide  range  of  values.  The  effects  of  these  variations  on  energy 
partitioning,  particle  shape,  depth  of  penetration,  and  target 
erosion  were  determined.  The  predictions  are  presented  by  graphs 
which  may  be  used  for  test  data  interpretation  and  extrapolation, 
test  planning,  hardware  design,  etc.  The  graphs  should  prove 
useful  because  one's  Intuitive  abilities  are  easily  overwhelmed  by 
the  complexity  of  the  impact  process. 
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